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Abstract 

We analyze how the range of disorder affects the localization properties of quasiparticles in a two- 
dimensional d-w&ve superconductor within the standard non-linear <r-model approach to disordered 
systems. We show that for purely long-range disorder, which only induces intra-node scattering 
processes, the approach is free from the ambiguities which often beset the disordered Dirac-fermion 
theories, and gives rise to a Wess-Zumino-Novikov-Witten action leading to vanishing density of 
states and finite conductivities. We also study the crossover induced by internode scattering due to 
a short range component of the disorder, thus providing a coherent non-linear cr-model description 
in agreement with all the various findings of different approaches. 

PACS numbers: 74.20.-z, 74.25.Fy, 71. 23. An, 72.15.Rn, 74.72-h 
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I. INTRODUCTION 



The low-temperature quasiparticle transport in two-dimensional c?-wave superconductors 
like cuprates is a fascinating issue due to the presence of four nodes in the energy spectrum 
of the Bogoliubov quasiparticles, around which the low-lying excitations have a Dirac-like 
dispersion. Within the self-consistent Coherent-Phase- Approximation in the limit of weak 
disorder, spin and thermal conductivities are found to be related by a Wiedemann-Franz 
law and to acquire universal values which do not depend on the disorder strength.- Inclusion 
of quantum interference in the framework of the standard non-linear a-model approach to 
disordered systems 2 -^ leads to a variety of universality classes^&^^^ii. In the "generic" 
case of short-ranged non magnetic impurities full localization of Bogoliubov quasiparticles 
is predicted- 

Nevertheless, with the only exceptions of YBCO(124)i2 and Pr2_ a: Ce a ;Cu04— , experi- 
ments in cuprates materials like YBCO(123}^^, BSCCO(2212) 15 and LSCO 17 do not 
show any evidence of strong or even weak localization in the superconducting phase down 
to 0.1 Kelvin. Various physical effects may be invoked to explain the disagreement between 
theory and experiments. 

For instance one may argue that the origin of the discrepancy are spin-flip scattering 
events, even though the systems are nominally free from magnetic impurities. Indeed, in 
the presence of spin-flips, the non-linear cx-model predicts that quantum interference has 
a delocalizing effect^. Alternatively, or in addition, strong dephasing processes might set 
the temperature scale for the onset of localization effects below the experimentally accessed 
region. This question has not been settled yet^. Residual interactions among quasiparticles 
can also favor delocalization&ii, even though, in the weak disorder limit, they are expected 
to be less effective since their coupling is proportional to the density of states which, already 
in the Born approximation, is very small. 

Another possible explanation invokes the range of the impurity potential. In the case of 
purely long range disorder, forward processes dominate and scattering occurs mainly within 
each node. In the extreme case of intra-node scattering only, it has been shown^22 that the 
density of states behaves quite differently from the isotropic-scattering case. In addition the 
eigenstates have been argued 20 to be delocalized, unlike for short-range impurity potential. 
Even though real disorder will always have an isotropic component which provides scattering 
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among all four nodes, and eventually drive the system to localization, one might argue that a 
large value of the intranode with respect to the internode scattering could lower the crossover 
temperature for the appearance of localization precursor effects. A sizeable amount of long 
range disorder has been indeed argued to be present in cuprates on the basis of STM and 
microwave conductivity experiments^. This is not surprising since superconducting cuprates 
are intrinsically disordered by the out-of-plane charge dopants which mainly provide a long- 
range disordered potential. Further doping with iso-valent impurities which substitute in- 
plane Cu-ions only adds a short-range component on top of the pre-existing long-range tail 
of the disordered potential. 

The results in the presence of purely intra-node impurity scattering have been 
obtainecr^i within an approach which is conceptually different from the standard non-linear 
cr-model approach to disordered systems. The latter starts from the Born approximation, 
i.e. from impurity-damped quasiparticles, and treats perturbatively what is beyond that, i.e. 
quantum interference effects on the diffusive motion. On the contrary, the alternative meth- 

rn 

ods used in Refs. [19|J20| do not rely on the Born approximation but just map the action of 
ballistic nodal-quasiparticles in the presence of disorder onto the action of one-dimensional 
(Id) fermions in the presence of an interaction, which is generated by the disorder average. 
Within the mapping, one of the two spatial dimensions transforms into the time coordinate 
of the Id model while the other into the Id spatial coordinate. The final model is then 
analyzed by abelian and non-abelian bosonization. The outcome of this analysis is that 
for purely intra-node or inter-opposite-node disorder, where essentially exact results can be 
obtainecU^^, the density of states is power-law vanishing at the chemical potential with an 
exponent which is disorder-dependent in the former case and universal in the latter. In both 
cases the results suggest that a diffusive quasiparticle motion never sets in, namely quasi- 
particles move ballistically down to zero energy. When the disorder also couples adjacent 
nodes, strong coupling arguments are invoked^ which predict localization of quasiparticles 
and linearly vanishing density of states. Yet, even this case seems to suggest a scenario in 
which quasiparticle motion from ballistic turns directly into localized without crossing any 
intermediate diffusive regime. 

The standard non-linear a-model approach applies the replica trick to average over dis- 
order from the outset. The resulting fermion interaction is then decoupled by introducing 
Q-matrix fields in terms of which an effective action is derived after integrating out the 
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fermions. The saddle point solution is just the Born approximation, which provides, in the 
case of nodal Bogoliubov quasiparticles, a finite density of states hence a finite damping. 
Finally, long-wavelength transverse fluctuations are taken with respect to the saddle point 
leading to a non-linear cr-model action for the Q-matrix fields. However, unlike in con- 
ventional disordered systems, in this particular case an additional term may appear in the 
non-linear cr-model action, namely a so-called Wezz-Zumino-Novikov-Witten (WZW) term. 
It was actually argued^^ that two opposite nodes share the same WZW term, while two 
adjacent ones have opposite WZW terms. As a consequence, when the two pairs of oppo- 
site nodes are uncoupled by disorder the WZW term is effective and the /^-function of the 
spin-conductance flows towards an intermediate-coupling fixed point, signaling a delocalized 
behavior. On the contrary, when disorder couples all nodes together, the two WZW terms 
cancel exactly and the non-linear cr-model has no more protection against flowing towards 
a zero-conductance strong coupling regime characterized by a linearly vanishing density of 
states^. 

From the above discussion one might be lead to conclude that the agreement between 
the conclusions drawn with either methods is merely an accident which does not justify per 
se the conventional non-linear cr-model approach when dealing with Dirac fermions. The 
main objection against the conventional non-linear cr-model is that it is not appropriate to 
start from a symmetry breaking saddle-point solution, associated to a mean-field-like finite 
density of states, when the outcome of including fluctuations is a vanishing density of states. 
Put in a different language, it is hard to believe in a method which starts by assuming a 
diffusive behavior if at the end it is discovered that a diffusive regime never appears. This 
criticism could invalidate also the results obtained with isotropic scattering, even though 
in this case the dimensionless coupling of the non-linear cr-model can be made small^ by 
assuming a large anisotropy of the Dirac dispersion (i.e. the velocities along and orthogonal 
to the Fermi surface). A related objection that can be raised is that the non-linear cr-model 
approach to disordered systems is commonly believed to be valid for length scales longer than 
the mean free path and energy scales smaller than the inverse relaxation time 1/r. On the 
contrary, the solution of the intra-node scattering problem demonstrates that disorder starts 
to affect for instance the density of states at energies of the order of the superconducting 
gap, hence much bigger than 1/r, namely in the regime when quasiparticle motion should 
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be still ballistic. 

It is the scope of the present work to clear up these inconsistencies of the non-linear a- 
model approach to disordered ci-wave supercondutors. We will demonstrate that the above, 
apparently contradicting, approaches can be actually reconciled. This is of particular interest 
since it provides further support to the standard non-linear cx-model approach based on the 
replica trick within the fermionic path-integral formalism, which remains so far one of the 
few available tools to deal with disorder in generic situations with a Fermi sea of interacting 
quasiparticles. Let us briefly summarize our main results. 

First we are going to present a simple and straightforward way to extract the WZW 
term. Indeed it is well known how to derive the WZW term within field theories defined on 
a continuous space with Dirac-like dispersing particles. However in disordered lattice models 
the existence of such a term is not at all a common situation. We will show that the WZW 
term emerges quite naturally within the conventional derivation of the non-linear cx-model 
for disordered systems as a consequence of the non-analytical properties of the spectrum 
within the Brillouin zone. More specifically, the spin-current density in momentum space, 
Jk, in a (i-wave superconductor is a 2 x 2 matrix in the Nambu spinor space. The WZW 
term arises just because the vector product Jk A Jk is finite and actually gives a measure 
of the vorticity around each node. In the presence of purely intra-node impurity scattering 
we obtain the same WZW action of the non-abelian bosonization from the Id mapping^>22, 
with however the inverse mean free path as a momentum upper cut-off instead of the inverse 
lattice spacing which is usual the Ultra- Violet-regularizer of the Id Dirac theory. In this 
context we elucidate the role of the WZW term in providing the correct scaling behavior of 
the density of states depending on the impurity-potential range. 

Another issue we clarify is the relationship between the coupling constant of the non- 
linear cx-model and the actual spin-conductance. The two quantities are known to coincide 
at the level of the Born approximation. However, rigorously speaking, the spin-conductance 
has to be calculated through a Kubo formula which involves advanced and retarded Green's 
functions. Since it makes a difference whether quasiparticles are right at the nodal points, 
with zero density of states, or slightly away from them, in which case the density of states 
is finite, one might wander whether the two quantities, spin-conductance and non-linear 
cx-model coupling, remain equal even beyond the Born approximation, especially in the case 
of intra-node and inter-opposite-node scattering. We will show that this is actually the case. 
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Finally we will show that the range of the impurity potential crucially affects the energy 
scale at which localization precursor effects starts to appear. In particular we will explicitly 
show that, keeping fixed the inverse relaxation time within the Born approximation, f/r , 
and increasing the relative weight of the long-range with respect to the short-range compo- 
nents of the disorder, leads to strong reduction below 1/tq of the energy scale for the on-set 
of localization. This may provide a natural explanation to the partial lack of experimental 
evidences of quasiparticle localization in superconducting cuprates. 

The paper is organized as follows. In Section II we introduce the model which we re- 
formulate within the fermionic path-integral formalism using the replica trick in Section III. 
The global symmetries of the action are discussed in Section IV. In Section V we start the 
derivation of the non-linear a-model which includes two terms. The "conventional" one is 
worked out in Section VI, and its drawbacks discussed in Section VII. The "unconventional" 
WZW term is derived in Section VIII and its consequences for intra-nodal disorder are 
discussed in Section IX. In Section X we analyse the role of inter-nodal scattering processes 
and finally Section XI is devoted to the concluding remarks. 



II. THE MODEL 



The model we study is described by an Hartree-Fock Hamiltonian for rf-wave supercon- 
ductors in the presence of disorder 

n = J>£ [ekT3 + A k 7i] ^ + J^)^T 3 ^ 
k kq 

= £4 hL 0) A + V ^ ^ r 3 A +g - (1) 

k kq 

Here V(q) is the impurity potential, tp^ = (cj^, c_ k |) the Nambu spinor, A k = A (cos k x a — 
cosk y a)/2 the superconducting gap with d-w&ve symmetry and e k the band-dispersion mea- 
sured with respect to the chemical potential. The Pauli matrices Tj's, i = 1,2,3, act on the 
Nambu spinor components. The spectrum of the Bogolubov quasiparticles is given by 

= y/<£ + A k , (2) 

and shows four nodal points at k x = k F (1,1)/V2, k- = — ki, k 2 = k F (— 1, l)/>/2 and 
k^ = — k 2 , k F being the Fermi momentum along the diagonals of the Brillouin zone. Around 
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the nodes it is actually more convenient to rotate the axes by 45 degrees through 

k^ = —j= (k x + k y ) , k v = -j= {k y — k x ) . (3) 

In the new reference frame, k x = kp (1, 0) and k 2 = kp (0, 1). In what follows we define 

* = k - k a (4) 

the momentum deviation from any of the four nodal points k a , a = 1,1, 2, 2. For small 
\k\ <C k F the spectrum around nodes 1(1) and 2(2) is, respectively, 



E kl+K = E 1 K ~ Jvp re| + v\ rf, (5) 



E k2+K = E 2k ~ yjvj? k\ + v\ «|, (6) 

thus having a conical Dirac-like form. Here vp = |VekJ and wa = |VAkJ. 

In the presence of disorder the motion of the gapless quasiparticles may become diffusive 
in the hydrodynamic regime. However diffusive propagators appear only in those channels 
which refer to conserved quantities. Hence, in superconductors, only thermal and spin 
density fluctuations might acquire a diffusive behavior. Let us therefore briefly discuss some 
properties of the spin current operator which are going to play an important role in our 
analysis. 

The ^-component of the spin-current operator in the Nambu representation satisfies 

Q ■ £ 4 Mq) A +q = £ 4 [< q - H®] A +q , (7) 

k k 

hence, for q — » 0, 

Mq) ^Jk = Ve k r 3 + VA k n . (8) 
Since the Pauli matrices anticommute, the following vector product turns out to be non-zero 

J k A J k = 2zr 2 Ve k A VA k . (9) 

The vector product (Q actually probes the vorticity of the spectrum in momentum space; 
nodes "1" and "1" have the same vorticity, «7 k A «7 k — > 2zt2VfVa, as well as nodes "2" and 
"2" , although opposite of node "1" , J k A J k — > — 2 % r 2 v F v&. We are going to show that this 
property is crucial to uncover the physical behavior in the presence of disorder. 
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III. PATH INTEGRAL FORMULATION 



To analyze the effect of disorder in this model, we are going to use a replica trick method 
within the fermionic path-integral approach. For that we associate to each fermionic operator 
Grassmann variables through 

C Lt ~~ *■ ^ko-? C k(J — > Ckcr- 

Notice that, unlike the original fermionic operators, Ck CT and Cko- are independent variables. 
After introducing the Nambu spinors 



^nk = (Ck|(^n), Qki (-«*>„)) , lp n]> 



where Ck a (iuj n ) and Ck a (iu} n ) are the Grassmann variables in Matsubara frequencies, the 
path-integral action without disorder reads 

So = ^2 t Ck 7-3 + A k n - iu n ) iPnk- (10) 

n k 

The disorder introduces the additional term 

Simp = V ( q ^ ^ r 3 T/Wq- (11) 

»t kq 

Since ip nk and are independent variables, the global transformation 

VVk^e^Vnk, t k -? nk e !flT2 , (12) 

becomes allowed within the path-integral and is indeed a symmetry transformation of the 
full action S = S + S imp when cj„ = 0. A finite frequency, u n 7^ 0, spoils this symmetry. If, 
in addition, the disorder is long range, namely it does not induce inter-node scatterings, it 
is possible to define four independent symmetry transformations of the above kind, one for 
each node. This type of chiral symmetry plays a crucial role for long range disorder, as was 
first emphasized in Ref. 

We notice that the disorder does not induce any mixing between different Matsubara 
frequencies, so that the action decouples into separate pieces, each one referring to a pair of 
opposite Matsubara frequencies, ±u n , which are coupled together by the superconducting 
term. For this reason we will just consider one of these pairs, with frequencies which we 
denote by ±uj, discarding all the others. This is enough to extract information about 
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the quasiparticle density of states (DOS) at given frequency as well as about transport 
coefficients. 

In order to derive a long-wavelength effective action, we resort to the replica-trick tech- 
nique, hence we first introduce N replicas of the Grassman variables, c^iza;) — > c a k,r(±iu;) 
and Ck a (±iu) — *> c ak(T (=ya;), with a = 1, . . . , N. At the end we shall send N — > 0. Next we 
define the column vectors Ck and Ck, with AN elements (N replicas, 2 spins and 2 frequencies) 
c k CT (±zo;) and c ako -(±za;), respectively. Finally we introduce new Nambu spinors through^ 

V2 y za 2 c k 

as well as 

* k = **_ k c t , 
with the charge conjugacy operator defined by^ 

c = ia 2 T\. 

Here the Pauli matrices <7i, g-i and 03 act on the spin components of the column vectors 
Ck and Ck, while n, r 2 and T3 act on the Nambu components (the diagonal elements refer 
to the particle- hole channels and the off-diagonal ones to the particle-particle channels). 
For later convenience, we also introduce the Pauli matrices s±, s 2 and S3, which act on 
the opposite frequency partners, as well as the identity matrices in the spin, Nambu and 
frequency subspaces, oq, r and sq, respectively. 

By means of the above definitions, the clean action can be written as 

S = t 6k +^ A kT 2 S! -iujs 3 ] tf k . (13) 

k 

Since we are interested in the low-energy long-wavelength behavior, we shall focus our at- 
tention only in small areas around each of the four nodal points. In the rotated reference 
frame the nodes lie at ki = k F (1, 0), k T = — ki, k 2 = kp (0, 1) and k 2 = — k 2 . Using the 
definition for the small momentum deviations away from the nodal points, we introduce 
new Grassmann variables defined around each node through 
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a — 1,1, 2, 2, as well as the corresponding Nambu spinors 



Co 



,*2« = I "'^ | , *2,t 



One notices that for a = 1, 2 




1 

i 



v c a,— ft) iCa,n cr 2) ^a— ft C ' 



The non-disordered action expanded around the nodes reads 

S = ^2 *i« + « «A S t 2 Si] - * I#C [uj? K f + iw A S r 2 Si] * Tfe 

+ *2k [«F Ki, + iVA K{ ?2 «l] *2ft - [ U F S + iv A T 2 Si] %fc 



(14) 



We find useful to define new spinors with components in each of the nodes through 

/*ift\ 



so that 



*2« 
V%ft/ 



*K V 1-ft' 2-ft' 1-tl U * -K n U ) 



where the Pauli matrices 7's act on the "m" and "m" subspace, namely connect two opposite 
nodes. This naturally leads to a new charge conjugacy defined through 



C = C71 = i<7 2 Ti7i. 



(15) 



In what follows we denote as pair 1 the two opposite nodes "1" and "1", and as pair 2 
the other two nodes, "2" and "2". Consequently we need to introduce also matrices in the 
two-pair subspace, 1 and 2, namely connecting adjacent nodes, which we will denote as po, 
the identity, p±, pi and p 3 , the three Pauli matrices. 
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The clean action reads now 

S = ^2 73 ^ + + "A72 si] ^ K 

+ ^nl3P3 ^ ( K ? _ S) faf - ^AT2«i] ^ K ~ ^ K S 3 $ K (16) 

K. 

= [H (K)-ius 3 ] V K (17) 

We notice that within this path-integral formulation, the four independent chiral symmetry 
transformations of the form (1121) translate into 

* -> Ttf, * ^ vj>CT*C*, 

with T given by any of the following expressions 

T = Q iSr2 Sl 70 p0 cr ' u (18) 

Ji — e *^-f2Si73Po 

y = gi S 7^ Sl 70 P3 0"-U (20) 

X 1 = e*^ r2 51 73 103 (21) 

where 9 is a phase factor and u an arbitrary unit vector. 

We assume that the scattering potential V(q) = V(— q)* induced by disorder has 
independent components which act inside each node, V(q) = Uo(q), between opposite 
nodes, V(q + 2ki) = U\{q) and V(q + 2k 2 ) = U 2 (q), and between adjacent nodes, 
V(q + ki — k 2 ) = Uu(q) and V(q + ki + k 2 ) = U l2 -(q), so that the impurity contribu- 
tion to the action is 



Simp = ^ U °( q ^ 

+ U M) ^Tn +q + Ui(q)* V jK+q 

K.,q 

+ ^ ^12(9) *2 K+q + % K %«+J + ^12(9)* \®2 K+q *1 K + %,c +q %J 

+ J] ^12(9) + *2k % K+9 ] + £/l2(9)* [% K+g *1k + *2*] • 

/t,qr 

(22) 
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Apart from the first term in the right hand side of ([22]) . we have for convenience kept the 
node labels. Since the following relations hold: 

UK. K 

E *2k%«+, = E = E *i«*w 

we can rewrite (122]) in the following way: 

S imp = E ^0(9) *K *»fq 

+E ^1(9) * 

+ E ^(9) ^ + ^2(9)* % K+ q ^2. 

+ 2 E ^12(9) *1« ^2 K +q + ^12(9)* ^2*+q *1* 

+2 E ^12(9) *i* %* +(J + ^2(9)* %* +(J (23) 

Finally we go back in real space, which now corresponds to the low-energy continuum 
limit of the original lattice model, and obtain the clean action 

S = J dr *(r) ^ 73 (-z0 € - ify) [v F + iv A r 2 si] *(r) 

+ y dr^(r)^7 3 p 3 (-id^ + id v ) [v F - i v A r 2 si] *(r) 

y dr W(r) s 3 \&(r) = y dr^(r) [H {r) - iu s 3 ) \&(r), (24) 

and the impurity term 

5 imp = y dr[/ (r)f(r)f(r) 

+C/ X (r) *i(r) %(r) + L7i(r)* %(r) ^(r) 

+f/ 2 (r) ¥ 2 (r) %(r) + U 2 (r)* %(r) * 2 (r) 

+2 C/ 12 (r) *i(r) * 2 (r) + 2 f/ 12 (r)* tt 2 (r) * a (r) 

+2£%(r) *i(r) %(r) + 2E/ l5 (r)* %(r) ^(r). (25) 
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Pauli matrices 


subspace of action 


CQ) &1, &2i &3 


spin components 




Nambu components 


SO, Si, S 2 , S3 


opposite frequency partners 


7o,7i)72,73 


opposite nodal points 


P0,Pl,P2,P3 


the two pairs of opposite nodes 



TABLE I: The various two-components subspaces in which each set of Pauli matrices act 
We further assume that the disorder is 5-like correlated with 



(Ua(r) 


> = 0, (U (r) 


Uo(r')) = 


u 2 5 (r 


-r'), 


(Ui(r) 


> = 0, (^(r) 


Uiir'T) = 


-- 2v 2 5 


[r - r') , 


(U 3 (r) 


> = 0, (U 2 (r) 


U 2 {r'Y) = 


-- 2v 2 5 


[r - r') , 


(Ui2(r] 


> = 0, (U 12 (r 


)MO*> 


= 2w 2 


5 (r — r' 




) = 0, (U 12 (r 


) U l2 {r'Y) 


= 2w 2 


5 (r — r' 



We conclude by noticing that (i) if only u ^ all transformations (fT8|) - (l2Tl) leave the 
action invariant; (ii) if u ^ and v ^ only (ITS]) and (120)) are allowed; (Hi) finally if u ^ 0, 
v 7^ and w^O only (ITS]) remains. 

IV. SYMMETRY PROPERTIES 

Since we have been obliged to introduce so many Pauli matrices, including the identity 
denoted as the zeroth Pauli matrix, for sake of clarity we prefer to start this Section by first 
summarizing in Tabled the subspaces in which any of them act. 

Now, let us uncover all global symmetry transformations. We start assuming that only 
intra-node disorder is present, namely U ^ but U\ = U 2 = Uyi = U 12 - = 0. 

If the frequency is zero, the action Sq + Si mp is invariant under unitary global transfor- 
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mations T such that 



They imply that 



as well as that 



We parametrize T as 



C'T'C^T = 73, 

C* T* C 73 r 2 si T = 7 3 r 2 Si, 

C t T t C'j3p 3 T = 73 p 3 , 

C t T t Cj 3 p 3 T 2 s 1 T = 73P3T2S1, 

C t T t C = T~ 1 ) (27) 
[T, 73 ] = 0, [T,p 3 ]=0, [T,r 2Sl ]=0. (28) 



T= l -(T 1+ T 2 ) p + l - {T x - T 2 ) p 3 , (29) 



where the suffix "1" refers to the pair 1 (opposite nodes 1 and 1) and "2" to the pair 2 
(opposite nodes 2 and 2). The symmetry modes for pair 1 and 2 can be parametrized by 
the 16iV x 16iV unitary operators 

1 1 / 14(2) 

Tl( 2 ) = - (70 + 73) 14(2) + 7> (70 - 73) 14(2) = 

2 2 v v i® 

where the V^'s are 8N x 8N unitary matrices. Because of ( 1271) we get 

nt rpt n rp—\ rp\ 

U 1 l(2) U — 1 l{2) ~ J l(2)> 

and 14(2) and 14(2) are actually not independent as 
Therefore 

Ti {2) = 



*i(2) 

c*(vj 2) )*c 

is indeed parametrized by a single 8N x 8N unitary matrix 14(2)- Since the two pairs of 
nodes behave similarly, in what follows we drop the suffices 1 and 2. According to ( 1281) we 
still need to impose that 

[V,r 2Sl ] = 0. (31) 
14 



To this purpose let us introduce the following unitary operation 



U 

which transforms 



1 , s 1 / 



e"^ S2 . (32) 



V->UV, V ^> V* U* c* = V* <* U* = V U* , 
so that for any operator 

One readily shows that 

[/t T2Sl [/ = S3) £^s 3 ?/ = -s 1 , (33) 

hence (EH) transforms into 



[V,s 3 } = 0, (34) 
which is fulfilled by the general expression 

V= l -(A + B) s ,+ l -(A-B) s 3 , (35) 

with A and B being independent 4iV x AN unitary matrices. Therefore the original sym- 
metry turns out to be U(AN) xU(4iV) for pair 1, and analogously for pair 2, in total 
U(4iV) xU(4iV) xU(4iV) xU(4JV). 

We notice that, in the presence of a finite frequency, uj ^ 0, we shall further impose 
through (1331) that 

[V,s 1 ]=0, 

which is satisfied by 

V = Cs , 

with C belonging to XJ(AN). Therefore the symmetry is lowered by the frequency down to 
U(4iV) for the pair 1 times U(AN) for the pair 2. 

Within the non-linear cr-model approach to disordered systems, the frequency plays actu- 
ally the role of a symmetry-breaking field. This leads to the identification of the transverse 
modes V 1 - as those satisfying 

Sx W- ax = (v a )- 1 = (v J -) t , 
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implying A = in Eq. (1351) so that we can write 



i 

So+ 2 



G 



S3, 



(36) 



where G belongs to U(4iV). In other words the coset space spanned by the transverse modes 
is still a group, namely U(4iV). It is convenient to factorize out of G the abelian component. 
For that we write 

G = eVf<* g, (37) 

where is a scalar and g belongs to SU(4iV). The forms (136]) and (1371) will be useful in 
the following to express the non-linear cr-model directly in terms of g and <fi. After the 
transformation (1321) from Eq. (I30|) we find that 



t 1 < 

n o 2 (V^ 2) ) n a 2 = [n o 2 si Vjjfi s x n a 2 ] 



which leads to 



; 1(2) 



) i(2), 9j(2) = n °2 [g 1(2) I 2 T X 



(38) 



(39) 



Let us conclude by showing what would change for more general disorder potential. First 
let us consider the case in which the disorder also contains terms which couple opposite 
nodes, i.e. U\ ^ and U 2 7^ 0. In this case the 73-modes are not anymore allowed by 
symmetry, so we have to impose that V\( 2 ) = Vj^, namely, through (1391) . that 0i(2) = and 

#i(2) = 0i(2) = r i a 2 g{ {2 ) °2 n ■ (40) 

The coset becomes now the group Sp(2iV) of unitary-symplectic matrices (also called 
USp(4iV)) for pair 1 and analogously for pair 2, i.e. Sp(2iV) xSp(2iV). 

Finally, if also scattering between adjacent nodes is allowed, Ui 2 7^ and U l2 7^ 0, then 
also the modes do not leave the action invariant. In this case we have to impose that 
Vf 1 = V 2 , so that the coset space is the group Sp(2iV). 

V. THE NON-LINEAR a-MODEL 

In what follows we begin analyzing just the case in which the disorder only induces intra- 
node scattering processes. At the end we will return back to the most general case. Within 
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S - U - 

t-'vmp 2 



the replica trick, we can average the action over the disorder probability distribution, after 
which (125]) with v — w — [see Eq. ff26l) ] transforms into 

S imp = -y j dr [nrMr)} 2 . (41) 
We define 32N x 32N matrices X(r) by 

X(r) = *(r)¥(r), 

so that (jJTJ) can be also written as 

J drTr(X(r)X{r)). (42) 

By means of an Hubbard-Stratonovich transformation one can show that 

exp j-y J drTr(X(r)X(r))j = J VQ(r) exp j- ^ rfr ^Tr [Q(rf] - z Tr [Q(r) X(r) 

(43) 

with Q(r) being hermitian 32N x 32N auxiliary matrix fields. 
In conclusion the full action, (I24"j) plus (143]) . becomes 

5 = ^/^Tr[Q(r) 2 ] 

+ / dr^(r) [H {r) -iQ(r) - ilu s 3 ] *(r). (44) 

One obtains the effective action which describes the auxiliary field Q(r) by integrating out 
the Grassmann variables, thus getting 

We further proceed in the derivation of a long- wavelength effective action for Q(r) by 
following the conventional approach. First of all we calculate the saddle point expression of 
Q(r), which we denote by Qos 3 , assuming it is uniform, in the presence of an infinitesimal 
symmetry breaking term 00 s 3 . Then we neglect longitudinal fluctuations and parametrize 
the actual Q(r) by 

Q(r) ~ T(r)- 1 Q s 3 T(r) = ~ (p + p 3 ) Qi(r) + ~ (p - Pa) Q 2 (r), (46) 
where through ([29]) Q a b(r) = 5 ab Q a {r) with 

Q a {r) = T Q (r)- x Q s 3 T a (r) = Q s 3 T a (r) 2 = T a ( r y 2 Q s 3 , (47) 
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being the auxiliary field in the pair subspace a, b = 1, 2. Notice that, even though the most 
general Q(r) would couple the nodes together, namely would include off-diagonal elements 
Qi2{ r ), (146 p does not contain any mixing term. This simply reflects that the off-diagonal 
components are not diffusive, hence massive. 

In order to derive the long-wavelength action we find it convenient to decompose the 
action into the real part 



S 



NLaM 



±- 2 J drTr[Q{rT] 



— -Tr In [iQ + iu s 3 — TCo) — -Tr In 



-iQ — iu s 3 — 7"4 



(48) 



which, as we shall see, gives the conventional non-linear cx-model, and the imaginary part 



Sr = — -Tr In [iQ + iu s 3 — 7Y ] + ^Tr In — iQ — iu s 3 — H 
which we will show gives rise to a WZW term. 



(49) 



VI. CONVENTIONAL cr-MODEL 



By means of fl4b|) . the real part of the action, fl4g|) . can be written as 

S NLaM = ^ 32iV Ql - ^Tr In [Q~ l ] , (50) 

where V e ff is the effective volume corresponding to the long- wavelength theory (roughly 
speaking V e // = V/4, since we have implicitly folded the Brillouin zone into a single quad- 
rant, in order to make all nodes coincide), and 

g- 1 = n Hi + Q 2 + uj 2 + uj{Qs 3 + s 3 Q)+iH Q-iQHi 

= Go 1 - X, (51) 

where 

g, 1 = n nl + Q 2 + u 2 , 

S = -u{Q,s 3 } -m Q + iQUl 
We notice that s 3 7i S3 = Hq, therefore 

i Ho Q - i Q H\ = iHo Q s 3 s 3 - i Q s 3 H s 3 = i [H , Q s 3 ] s 3 = J ■ VQ, 
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where we made use of the equivalence Q s 3 = T 2 and of the expression of the spin current 
operator in the long wave-length limit 



J(r)=i [H {r),r]. 



(52) 



Analogously one can show that 

iHoQ — iQW — is 3 Hls 3 Q — is 3 s 3 QTil = -is 3 s 3 Q,Hl = VQ ■ J\ 
so that the self-energy operator can be written as 

£ = —to {Q, s 3 } - J ■ VQ = -to {Q, s 3 } - VQ • J f . 



(53) 



We then expand the action up to second order in E and obtain 
S NLaM = ^32iVg2_I T r ln^- 1 ] 

= ^ 32iV Ql + ^Tr In [Go] - ^Tr In [1 - G E] 

= ^ 32iV g 2 + I T r In [0 O ] + ^Tr [Q E] + ^Tr [0 O E £ S] 

= ^§ 32iV g 2 + ^Tr In [&] - |Tr [Q s 3 So] 

+ ^Tr [0 O {g, s 3 } Go {Q, s 3 }] + ^Tr \g VQ-J^GoJ- VQ 



(54) 



Notice that we have arrived to SnlvM in terms of VQ without passing through a gauge 
transformation to carry out the gradient expansion so avoiding any problem related to a 
proper treatment of the Jacobian determinant 2 - . 



A. Saddle point equation 

In momentum space one finds that 

fto«?4 K = \ { K l + K l) [ V F + ^a) + P3 7j (4 - K%) (v 2 F - vl) 

= \ (K. + El K ) + P3 \ (El - E^) , 

where E ln is the spectrum of quasiparticles around nodes "1" and "1" , and analogously for 
E2K, see (|5]) and ([6]). Therefore 

Gok, = 77 [Gin + G2k] + P3 77 [Gin — Gln\ , 
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where 

a 1 a 1 

EL + Ql + ^ y2K EL + 

so that 

GL = \ [GL + G L] + Pa \ [GL - GL] (55) 

The saddle-point equation is obtained by taking Q(r) = Q and minimizing the action, 
which leads to the self-consistency equation 



-^E F 2^^ ,2 > (56) 



Q + u ■ V ^ El + (Q + uj) 
originally derived in Ref. The solution of this equation, in the limit uj — > 0, reads 

/ 7TV F V A \ 

exp (~2^-J , 

where A is an ultraviolet cut-off which is roughly the energy scale above which the spectrum 
deviates appreciably from a linear one, in other words A ~ A. We notice that in the generic 
case with finite inter-nodal scattering, namely with non zero v 2 and w 2 , the saddle point 
equation remains the same apart form u 2 — > u 2 + 2w 2 + v 2 . 

The quasiparticle density of states N(e), after the analytic continuation of iuo to the 
positive real axis, iu — >• e > 0, turns out to be, for small e, 



JV(c) 



7T 2 VfVa 



Qo In . + tan" 1 



e 2 + Q2 2 V2 2eg 



(57) 



Therefore the density of states at the saddle-point acquires a finite value ~ Qo/u 2 at the 
chemical potential e = 0, while, for e ^> Qo, turns back to the linear dependence N(e) ~ e 
as in the absence of disorder. 

B. Frequency dependent terms 

The first order term in the frequency is just 

- |Tr [Q s 3 Go] = ^ J dv Tr [Q(r) s 3 ) . (58) 
The second order term is readily found to be 

yTr [G {Q,s 3 } G {Q,s 3 }) = / dxTr [{Q(r),s 3 } 2 ] . (59) 
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This term is negligible as compared to ([58]) for frequencies u much smaller than vpv^No- 
The latter is therefore the energy scale below which diffusion sets in at the mean field level, 
namely the so-called inverse relaxation time in the Born approximation, l/r . 



C. Gradient expansion 



The Fourier transform of the current operator is 

J{q) = 7s g WQ) + Mo)] + 7s Ps ~ [Ji{q) ~ Mq)) , 
where, at long wavelengths q — > 0, 

Ji(q) -> Ji = (v F , iv A r 2 s 1 ) 
J 2 {q) -> J 2 = {iv A T 2 sx , v F ) . 

It is then easy to derive the following expression of the tensor product J ® 
J®J ] = 



J i j * j \)=i(v F +vD i10 



1 



1/2 2X / 1 \ 01 

2 » -1 / -10 



To evaluate the second order gradient correction in (|54|) we notice that 
^Tr Vg • J 1 " £0 J • VQ] = i ^ Tr [$Q fyQ £ ■/] £° J * 
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^ Tr \diQ d 3 Q (g j] g Ji + g j\ g j 3 



(60) 



(61) 



l-.l 



In the long wavelength limit 

go Jj go Ji + go jJ Go Jj — y go*, j] go*. j% + go*. j\ go*, 
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since £ K g 2 ln - G 2 K = 0. We find that 



1 /" c/k;^ 



i r°° . 2 / i \ 2 i i 



lnv F v A Jo \x 2 + QIJ 8nv F v A Qo 

The Drude spin- conductivity is defined byi 



47T 2 fi?fA 

so that the second order gradient correction can be written as 



/2 2 
drTr[Vg(r).VQ(r)]+^^Tr{p 3 [^g(r)^g(r) - d v Q(r) d v Q(r)]} . 
v z F + v A 



D. The final form of Sjvlo-m 



Collecting all contributions which are second order in the gradient expansion and first 
order in the frequency we eventually obtain 



S NL aM = Y^Ql a J drTr[VQ(r) -VQ(r)} 



77 ° [ dr^—^Tr{p 3 [d^Q(r)d^Q(r) - d v Q(r) d v Q(r)}} 
J v F -+- v A 



ieg, 



■' ' /•■ " ' ' A 



-uJ^JdrTr[Q(r)s 3 ] 

i=l,2 ^° 



^^Tr[Q,(r) S3 ]. (63) 



Here we have defined a metric tensor for the pair 1 which includes nodes "1" and "1' 



cr 



2 I v 2 p 

V 2 F + V\ \ (J v l 



and for the pair 2, i.e. for nodes "2" and "2" 

2 / ^ 



a m - i \ 



^F + V A \ V 2 



F 
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In reality the above action is not the most general one allowed by symmetry As we 
discussed earlier, the theory possesses two chiral abelian sectors, which actually occur in the 
singlet channels 73T 2 Sip an d 73T2S1P3, see ( TH)|) and (12"T1) . In analogy with Ref. [241 ]. we 
expect that upon integrating out longitudinal fluctuations the following term would appear: 

5Snl*m = an V / drTrfrsTiSxQiir^Qiir)} a^ r Ii[^T 2 s 1 Q i {r)d v Q l {r)}, 

(64) 

with II oc u 2 . 

In conclusion the most general non-linear a model is given by 

S nl *m = Yl [dr-^aTrld.Q^a^dMr^-u^TrlQ^Ss] 

+ -^ami[ l3 r 2 s 1 Q t (r)d,Q l (r)] a® Tr[y z T 2 s 1 Q i (r)d„Q i (r)) . (65) 

VII. FAILURE OF THE CONVENTIONAL cr-MODEL 

The non-linear a-model ( l65l) belongs to one of the known chiral a-models encountered 
when two-sublattice symmetry holds, see Refs. jflfl. If we simply b orrow known results, 



see e.g. Table I in Ref. [9|, we should expect that the /3-function of the conductivity vanishes 
in the zero replica limit, N — > 0. That would imply absence of localization and persistence 
of diffusive modes. Moreover we should predict a quasiparticle density of states N(e) which 
diverges like^^^ 

1 r . / r 

(66) 



N(e) ~ - exp 
e 



-Ax 3 1 In- 



with A and B some positive constants. This is clearly suspicious since in the absence of 
disorder the density of states actually vanishes, N(e) ~ e. 

One is tempted to correlate the above suspicious result with what is found for the ele- 
mentary loops of the Wilson- Polyakov renormalization group approach. Here one integrates 
out iteratively degrees of freedom in a momentum shell from the highest cut-off A to A/s, 
with s > eventually sent to infinity. In our case these fundamental loops for either nodes 
are given by 

If dK 1 

= In s = g In s, 



TXO 



A/s<|k|<A a^u K v 
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and provide the definition of the dimensionless coupling constant g, which should necessarily 
be small to justify a loop-expansion in g In s <C 1. In our case it turns out that g — 1, making 
any loop expansion meaningless. 

This results is at odds with the standard non linear a-models for disordered systems 
where g ~ 1/cr <c 1 for weak disorder. Here, whatever weak is the disorder, yet g = 1. This 



peculiar fact was originally discussed in Ref. 19j, where the authors identified correctly the 



complete failure of the non-crossing approximation as a starting point of perturbation theory 
due to the absence of small control parameters. This might lead to the conclusion that the 
non-linear cx-model we have so far derived is useless in this problem, since it heavily relies 
upon the assumption that quantum fluctuations around the saddle-point solution can be 
controlled perturbatively. 

In the following Section we are going to show that the above conclusion is not correct and 
that one only needs to be more careful in deriving the proper non-linear a-model action. 

VIII. WESS-ZUMINO-NOVIKOV-WITTEN TERM 

Indeed we have not yet accomplished all our plan, as we still need to evaluate the imag- 
inary part of the action (1491) . At leading order we can drop the frequency dependence of 
fll9|) . hence 

S r = - jTr In [iQ - H ] + ^Tr In [-iQ - wj] . (67) 
It is more convenient to evaluate the variation of along a massless path defined through 

5Q{r) Q{r) + Q{r) 5Q{r) = 0. (68) 

We find that 

SS r = -'-Tr [G 8Q\ + ^Tr [G f SQ] , (69) 

where 

G=(iQ-H y 1 . (70) 
We notice that, in the long wavelength limit, 

(-iQ - 74) {iQ - H ) ~ — VQ ■ J — - J ] ■ VQ, (71) 

and 

{iQ - Ho) (~iQ - H\) ~ Go 1 + VQ • J = Go 1 + J ] • VQ, (72) 
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hence 

g = [i - g vq ■ J}- 1 Go (-zQ - ?tf 

~ [1 + Go VQ ■ J + Go VQ ■ jg VQ • J] Go [-iQ - Wj) , 
Gt = [1 + Go VQ ■ J]' 1 Go (iQ - H ) 

^ [1 - Go VQ ■ J + Go VQ ■ JGo VQ ■ J] Go (iQ - Ho) . 

The leading non- vanishing contribution to (1691) derives from the second order gradient cor- 
rection to G — G^, which reads 

G-G^ ~-2iG VQ- JGoVQ- JGoQ- (73) 

Therefore 

SS r = ~Tr [5Q Go VQ ■ J Go VQ ■ J Go Q] 

= -~Tr [SQ Go VQ ■ J Go J ] • VQ Go Q] (74) 

The only term which contributes comes from the anti-symmetric component of the tensor 
J®J ] : 



1 



- \ J { Jj - Jj Jj j = -i tij v F v A r 2 si p 3 , (75) 
being the Levi-Civita tensor, thus leading to 

5S r = % r %^ fijQ Q 6 Q d *Q T ^ s i ps] ( 76 ) 

where 



K 

OO 



If fly 11 

= 8W dx2 W + Ql) = i^Ql' (7?) 

Introducing a fictitious coordinate which parametrizes the massless path, we finally get 

Sr = 1 !<L Wo I d3r e " vv Tr [r2 si p3 g(r) ^ g(r) g(r) dpQ{r) g(r) dqQ{r)] (78) 
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This actually represents opposite WZW terms for each of the two pairs of nodes, namely for 
each of the two independent U(l) xSU(4iV) a-models. It is clear that if the disorder coupled 
the two pairs of nodes, we would be forced to identify Qi(r) with Q 2 (r), so that the WZW 
term would cancel in that case. To make (1751) more explicit we remind that, for a = 1, 2, 

Q a (r) = T a (r)- 2 s 3 Q = s 3 Q T a (r) 2 , (79) 

where 

T>) = ~ ( 7o + 73) KV) + ~ (7o - 73) c* [V7 2 (r)f c, (80) 



and 



~ (s + s 3 ) e J v/# Mr) ga{r) + ~(s - s 3 ) eVF'-W g a {r)^ 



U\ (81) 



KT(r) = U 

U being defined in Eq. fl32|) . By means of fl79|) . flHUl) and flHTj) . the expression flTHl) can be 
finally written as 

<S r = y rf 3 f e^Tr [g x (r) ] d^r) gi(r) j d^r) gi (r) ] d^r)] 

-i ^ J d'r e^Tr ^ 2 (r) t ^ 2 (r) g 2 {r)^ d u g 2 {r) g 2 (r)^ d v g 2 (r)} , (82) 

appropriate for two WZW models SU(4iV) 1 . 

Let us express all other terms in the action by means of <p a and g a . First of all the density 
of state operator becomes 

1 2 

— Tr [s 3 Q(r)] = Tr [T 2 (r)j = 2 £ e* V^^Tr [<?„(*•)] + e" 4 Vf ^Tr [^ a (r)t] . (83) 



a=l 



Qo 

Moreover one readily finds that 

Tr [d M Q (r) a$ d u Q a (r)] = 16n d^ a {r) otjg ^0 a (r) + 4Tr [9^ a (r) a$ 0„<7 a (r)t] 



and 



^Tr [73 T 2Sl Q a (r)d,Q a (r)} a$ Tr [ l3 r 2 Sl Q a (r)d u Q a (r)} = 



In conclusion the action expressed in terms of the fundamental fields and including the 
WZW term reads 
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Swzw = fdr^aTr [d,g a (r) d v g a {r) ] ] + n 2 a (1 + NU) d^r) a<$ d^r) 

a=l,2 J 

-lottN je^^WTr [g a (r)\ + e'^^Tr [g a (r)*] } + 4°\ (84) 

with 

4 a) = ±i ^i- J d 3 r e^ v Tr [g a (r)i 9^ a (r) </ a (r)t d„<7„(r) </„(»•)* ^ a (r)] , (85) 
where the plus refers to a = 1 and the minus to a = 2. 

IX. CONSEQUENCES OF THE WZW TERM 

We have just shown that the actual field theory which describes <i-wave superconductors 
in the presence of a disorder which only permits intra-node scattering processes is not a 
conventional non linear cr-model but instead it represents two decoupled U(l)xSU(4JV)i 
WZW models. Moreover if, for instance, we consider pair 1, then 



7T 2 

4 

7T 1 V 2 F + V 2 A 2 



a HrTHr I d " r ^Tr [% a (r) % a (r)t] + v 2 A Tr [d v g a (r) 9^ a (r)t] 

V F ' V A J 

d 2 r v 2 F Tr [% a (r) % a (r)t] + v 2 A Tr [d v g a (r) d v g a (r)^ 



4 47T 2 VfV^ v 2 f + v A 

which, upon the change of variable £ — > i>.f£ and 77 — > i>a?7, shows that each WZW model is 
right at its fixed point. Hence there is no ambiguity in the zero replica limit. 

Now we can draw some consequences of what we have found. The first is that the average 
value of the density of states (DOS) at the chemical potential stays zero, as in the absence of 
disorder and contrary to the Born approximation. In particular the dimension of the density 
of states operator in the zero replica limit iV — > is 

N-l 1 

a = + — -> i - n, 

while the dimension of the frequency is [uo] — 2 — (1 — IT) = 1 + II. This implies that the 
DOS at finite frequency behaves as 

N(u) ~ Mmr, (86) 
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in agreement with Ref. 19J. 



Notice the fl86l) stems from the fact that the WZW terms modifies the non-linear cr-model 
results leading to Eq. (1661) in two ways: i)it makes II unrenormalized as well as cr, ii) it 
adds the further contribution one to the dimension of the density of states operator. The 
only difference between the Id mapping and the non-linear cr-model results for the density 
of states are the constant factors fixed by the range of validity. Id mapping: N(u) ~ 
l/(v F v A )iu(iu/A)- m ^ 1+ ^ with uj < A; non-linear cr-model: N(uj) ~ iV (Wo) (1 - n)/(1+n) 
with uj < l/r , iV being the saddle point value of the density of states. At leading ording 
in u 2 the matching of the two expressions is provided by NqTq ~ 1. 

The second consequence concerns the transport properties. We have shown that within 
the Born approximation a finite spin-conductivity arises. Is this result still true beyond that 
approximation? The renormalization group analysis which identifies the spin conductivity 
with the coupling of the non-linear cr-model would say that a s stays unrenormalized to its 
fixed point value. Can we understand this in the Id mapping language?. 

In order to answer this question, we introduce a uniform spin vector potential within the 
action, which in our path-integral approach has the form 

A = A cr 3 + A 1 a 3 s 1 . (87) 

The action becomes now a functional of A, i.e. S — > S(A), and the spin-conductivity turns 
out to be given by 




dA lv ' A. I) 



where 



Z{A) = J DQ(r) 



is the generating functional in the presence of A. It is possible to show that at second order 
in A one needs just to make the following substitution in the action 

^Q(r)->0 M Q(r) + I[Q(rM M ]. (89) 

Upon the action of fl32|) 

A -> U ] A U = A <7 3 + Ai cr 3 T 2 s 3 , 



28 



and we can show that (189]) implies the following transformation law of the matrix field G(r) 
in the presence of A 

G(r) — > 6^5 / ' dr'-{A (r')a z -A 1 (r')T 2 cr 3 ) Q( r } e i I dr' ■ (A (r') a :i +A 1 (r') r 2 a 3 ) ^ ^qq^ 

In order to better understand the role of A it is convenient to translate the WZW action 
into the language of the underneath free one-dimensional fermions. One may identify the 
component Gi ;a p(r) of the matrix field for pair "1", where a and (3 run over AN indices, 
with 

Gw(^i^ a W*wW> (91) 

being ^i- a R(r) and $i ;a i(r) right and left moving one dimensional Fermi fields, respectively, 
and the two component vector r playing the role of space and time coordinates. Since pair 
"2" has the opposite WZW term of pair "1", it is more appropriate to define 

G 2 Mr) ^ -i%. La (r)*i >R/3 (r), (92) 



which formally yields to equal WZW terms. Then ( 1901 implies for the fermions the trans- 
formation laws 

Jr) — > e~^ fdr'-(Ao(r')a 3 +Ai(r')v2CT3) (7-^ 

^ 2 . R (r) ^ e-^ dr '- {Ao(r ' )a3+Al(r ' )T2a3) * 2;fl (r), 
which we may interpret as if A couples to the spin-density operator 

i=l,2 

while A\ to the spin-current operator 

*U( r ) a3 T2 *i;*( r ) ~ *U( r ) a3 r2 *i;i( r ) " ^U( r ) ^ r 2 ^ 2;i? (r) + ^ 2;L (r) a 3 r 2 * 2;£ (r). 

(94) 

Since the fermions are free, apart from the abelian sector which is not coupled to A, the 
susceptibility towards Aq is the same as that towards Ai, hence the spin-conductivity would 
seem to vanish at u strictly equal to zero, again unlike what we found within the Born 
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approximation. Actually one has to be more careful in drawing such a conclusion. Let us 
suppose to do the same calculation at finite frequency uj and afterwards send uj — > 0. A 
frequency uj 7^ plays the role of an explicit dimerization in the one dimensional fermionic 
problem: 

-ttuNo Tr [G.(r) + Gj(r)] = -*™iV £ *U(r)* i;I , a (r) - *1 ;L »*. fla (r). 

i=l,2 i,a 

(95) 

It is straightforward to show that in the presence of a finite dimerization the current-current 
susceptibility is finite and practically equal to that one in the absence of dimerization, 
while the density- density susceptibility is zero. This discontinuous behavior at uj = as 
opposed to uj —>■ is again a manifestation of the chiral anomaly which plays such an 
important role in this problem^. Since it is physically more appropriate to identify the 
spin conductance through the uj — > limit, we conclude that, in spite of the vanishing 
DOS, the spin conductivity is finite. In other words, in spite of the fact that the DOS is 
vanishing at the chemical potential, namely that quasiparticle motion is undamped hence 
remains ballistic, yet the spin conductivity acquires a finite value in agreement with the 
Drude approximation. 



X. INTER-NODE SCATTERING PROCESSES 

So far we have just considered the role of impurity scattering within each node. Now let 
us extend our analysis by including also inter-node scattering processes. Upon integrating 
out the most general disorder, we find two additional terms. The first describes opposite 
node scattering processes, and it reads: 

SSL P = -v 2 J dr p?!(r) %(r)] [%(r) ^(r)] + [* 2 (r) %(r)] [%(r) tf a (r)] . (96) 

The second derives from the impurity scattering processes which couple the two pairs of 
nodes, and it is given by 

«C = -w 2 J dr P?!(r)* 2 (r)] [* a (r) *i(r)] + [%(r) %(r)] [% (r) * T (r)] 

+ [* x (r) %(r)] [f a(r) ^(r)] + * 2 (r)] [* 2 (r) %(r)] . (97) 

As before we can decouple the four-fermion terms by introducing auxiliary Hubbard- 
Stratonovich fields. Since these fields are expected to be massive, we can further expand the 
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action up to second order in those fields after integrating out fermions. We finally integrate 
on the auxiliary massive fields. The net result is still an action for the Q matrices only 
which includes now the additional terms: 

N f 

6SL P = - X 1 ^ J dr Tr [Q(r) 7l Q(r) 7l ] , (98) 



and 



where 



0C> imp — ~ A 77" 



f dr ^ Tr [Q(r) Pl 7 Q(r) Pl 7< ] , (99) 

i=0,l 



X 7 ^ x/1 ^ 
A ~ — , A ~ — . 



The first term ([98]) tends to lock G x = G- and G 2 = G^, while (JHU) locks Gj. = G 2 . 

When opposite node scattering is added, still leaving pairs of opposite nodes uncoupled, 
only the symmetric Q-combinations of opposite nodes stay massless, the Il-term disappears 
and the coset for soft modes is Sp(2N) for each pair of nodes. The /3-function is vanishing 
only because of the contribution of the WZW term, and density of states vanishes with a 
universal exponent in agreement with the known result In the absence of isotropic 
scattering the vanishing of the (3 function still indicates that the spin and heat conductivities 
have a metallic behavior. 

Finally, in the generic case in which all nodes are coupled only the four nodes symmetric 
combination of the Q's is required for the soft modes. The coset is again Sp(2N), but now 
it represents degrees of freedom coming from all nodes. The two WZW terms are written in 
terms of the same Q and they cancels since they have opposite sign. The action then reduce 
to the Snl<tM as derived by [5(. 



A. Scaling analysis of the general model 

In order to elucidate the role of inter-node scattering processes, it is convenient to trans- 
form the 2 dimensional non-linear cx-model into a 1+1 dimensional model of interacting 
fermions. 
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a a 



i,x,a, a i,x,G, a l,x,a,a %X,G, a 

" X " 

i ,t ',o ', a' i ,x ',a a 2,x ',a ', a ^ ^ l,x \o >, a' 

§ g 3b 

i,x,a,a i,x',a',a' l,x,a,a l,x',a',a' 

i,-x,-a,a ^i-x',-a',a' G<7 2,-x,-a,a J 2,-T',-CJ',a' 

- - -► - - - > fc > 



FIG. 1: Interaction vertices for the most general disorder. Solid(dashed) lines refer to right(left) 
moving fermions. The label % = 1, 2 refers to the two pairs of opposite nodes, r, r' = ±1 to the two 
Nambu components, a, a' =f , j to the spin and a, a' = 1, . . . , N to the replicas. The symbol <r <r' 
in front of the and 535 coupling constants means +1 if a = a' and -1 otherwise. 

We represent the matrix fields for pair "1" according to: 

Gl; TC ra,T> a> a> = ' '^1;R,t a a ^1;L,t' <t' a' ' 
^TiTiijT-'ir'a' = — * ^l;L,raa ^ 1;R,t' a' a' > 

G i ;WlT v. = (^i^Gl^n) (100) 

' ' \ / r a' a, r a a 

1;t a a , r a' a' L,n,—T a a t a a 

As we mentioned, since pair "2" has the opposite WZW term of pair "1", it is convenient 
to define G 2 . Taa T , a , a , = —i^^-L Taa^l-Rr' a' a'- m sucn a wa y the two WZW terms become 
equal leaving no ambiguity in mapping the non-linear a-model onto a one-dimensional model 
of interacting electrons with the interaction vertices drawn in Fig. [TJ We notice that the 
coupling g in Fig. [1] derives from the two terms in Eq. flM|) . The bare values of the coupling 
constants are approximately 

(0) (0) 2 
9 b = 9 3b - w ■ 

Generally u 2 > w 2 > v 2 , the equality holding only for short range impurity potential. 
An important observation is that two-loop corrections to the renormalization group (RG) 
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equations vanish in the zero replica limit, hence the RG equations valid up to two-loops are 
found to be: 

d9 2,2,2 

9s+9b+ 9zb, 



dins 

dg 3 

dins 
dg b 

dins 
dgzb 

dins 



4=993 + 4g b g 3b , 
2ggb + 2g 3 g 3b , 
2 £3 9b + 2g 3b g, 



where s — > oo is the scaling parameter. As discussed in Ref. [20J, the velocity anisotropy 
does not enter the RG equations, which remains true at least up to two loops in our fermionic 
replica trick approach. It is convenient to define g± = g 3b ± g b , so that 

dg 2 i 1 / 2 i 2 \ 

9 3 + o \9+ + 9-) , 



dins " 2 

dg 3 



dins 

dg+ 

dins 

dg- 

dlns 



4gg 3 + (g 2 + - gl) 
2 (g + g 3 ) g+, 
2 (g - gs) g- 



Given the appropriate bare values of the amplitudes, one readily recognizes that the RG 
flow maintains the initial condition g_ = 0, hence the scaling equations reduce to 

-9i + l9l, (101) 



dins ^ 2 

^ = ±9 93 + 91 (102) 
dins 

= 2(g + g 3 )g + . (103) 



dg+ 



dins 

The RG equations with the appropriate initial conditions always flow to strong coupling 
with 

1 1 

93 ~ % 9 ~ 9+ -> r i 1 — , 

3 In s c — In s 

where s c can be interpreted as the correlation length of the modes which acquire a mass gap 
E c by the interaction, with 

E c * (104) 



Therefore the two pairs of nodes get strongly coupled, in agreement with Refs. 19lJ2Q|]22 ]. 
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B. Strong coupling analysis 



In order to gain further insight into the strong coupling phase towards which RG flows, 
let us consider the case of a single replica N = 1. For further simplification it is convenient 
to adopt the same approach as in Refs. 019 20 22 and neglect the role of the opposite 



frequencies, which amounts to drop the r-label. The model thus reduces to two interacting 
chains of spinful fermions, each chain representing a pair of nodes. The coupling g of Fig. [1] 
only couples to the charge sector and makes the intra-chain umklapp, the coupling g%, a 
relevant perturbation which opens a charge gap on each chain. Therefore the model is 
equivalent to two coupled spin-1/2 chains. If we denote by n\^{x) and ex( 2 )(x) respectively 
the staggered magnetization and the dimerization of chain 1(2), the coupling among the 
chains is ferromagnetic and given by 

w 2 / dxei(x) e 2 (x) — Ani(x) ■ n 2 (x). 



As shown in Ref. |29_|, this model is equivalent to an SO (4) Gross- Neveau which turns out to 
be fully massive or, equivalently, by four two-dimensional off-critical classical Ising models, 
three ordered and one disordered, (01) = (cr 2 ) = (er 3 ) = (/i 4 ) = <j / 0, where Oi and /ij, 
% = 1, ... ,4, are order and disorder parameters, respectively. The ground state is rigid to 
an external magnetic field and to a spin vector potential opposite for the two chains, hence 
the conductivity is zero. As we discussed, a finite frequency amounts to add a term 

lu / dxei(x) + e 2 (x) oc uj I dx o~\(x) cr 2 (x) as(x) cr 4 (x) ~ u cr 3 / dxa^x), 



which actually plays the role of an external magnetic field acting on the fourth disordered 
Ising copy. The net result is that (ct^x))^^ ~ to, which in turns mean that the DOS remains 
linear in frequency. Even in the presence of an explicit dimerization, the susceptibilities 
towards a magnetic field or towards a spin vector potential opposite for the two chains still 
vanish. Were these results valid for any N, even for iV — > 0, we should conclude that i) the 
model is indeed insulating; ii) the DOS is linear in frequency, in full agreement with Ref. 



.B 



C. Identification of relevant energy scales 



The previous analysis of the N = 1 model shows that the vanishing of thermal and spin 
conductivities in the model for a disordered d-wave superconductor translates in the language 

34 



of the effective one-dimensional fermionic model into the existence of a finite spin-gap. The 
correlation length associated with this gap should then represent the localization length of 
the Bogoliubov quasiparticles. We may estimate this correlation length as the scale s c at 
which the RG equations (I101H1031) encounter a singularity. In addition we may introduce 
the mass gap E c through Eq. (11041) which can be identified as the energy scale around which 
localization effects appear. It is worth noticing that s c underestimates the spin correlation 
length, hence the actual localization length. The reason is that the RG equations blow up 
on a scale which is related to the largest gap in the excitation spectrum. Since the coupling 
g only affects charge degrees of freedom, the largest gap is expected in the charge sector, 
the spin gap being smaller. Keeping this in mind, in what follows we shall discuss how s c , 
or better E c , depend on the range of the disorder potential. 

First of all we need to identify some reference scale to compare with E c . The natural 
candidate would be the inverse relaxation time l/r in the Born approximation. In the 
generic case u 2 >w 2 > v 2 , l/r = 2Q where Qo is obtained by Eq. (l56l) with u 2 substituted 
by u 2 + 2w 2 + v 2 . We expect that the actual E c is always smaller then 1/ro, the two values 
being closest for extremely short-range disorder, namely u 2 = w 2 = v 2 . Since the derivation 
of the non-linear a-model does not provide with the precise dependence of the initial values 
of the coupling constants, g(°\ g% and g+ , on the impurity potential, we will assume that 
the short range disorder corresponds to 2g^ = gf^ = g+ and moreover that the mass gap 
E c in this case can be identified as 1/ro i n the Born approximation. Upon integrating the 
RG equations (1101111031) with this initial condition, the value of E c ~ l/r is found to be 

^M^^f+j' <105) 
where A is the ultraviolet cut-off of the order of the gap A. In order to appreciate the role 
of the range of the disorder-potential, let us analyze the RG equations (I101H103I) keeping 
fixed the combination g = 2g^ + g^ + g+\ i.e. at constant 1/t , and increasing the value 
of g(°) at expenses of gf^ + g^P = go — 2g(°\ One readily realizes that as g^ increases 
E c decreases from its short-range value (I105p . For instance, if we assume gfg = and 

gf = 9o- 2 s (0) < 9o. th ™ 

To y 4g J T ' 

which explicitly shows that localization effects may show up at energies/temperatures much 
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smaller than the 1/tq in the Born approximation. 



XI. CONCLUSIONS 

In this work we have presented a derivation of the non-linear a-model for disordered 
d-w&ve superconductors able to deal with a finite range impurity potential, namely with an 
intra-node scattering potential generically different from the inter-nodal one. Within this 
derivation we have been able to clarify some controversial issues concerning the validity of 
a conventional non-linear a-model approach when dealing with disordered Dirac fermions. 
We have indeed found that the non-linear q-model approach is actually equivalent to the 
alternative method first introduced in Ref. 
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which consists in mapping the disordered 



model onto a one-dimensional model of interacting fermions. The energy upper cut-off 
is provided by the inverse Born relaxation time l/r in the non-linear cr-model approach 
and by the superconducting gap A in the Id mapping. A closely related aspect which also 
emerges for d-wave superconductors is the existence of a Wess-Zumino-Novikov-Witten 
term related to the vorticity of the spectrum in momentum space. Both these features are 
responsible of several interesting phenomena. For instance, unlike conventional disordered 
systems, in this case disorder starts playing a role (particularly in the density of states) 
when quasiparticle motion is still ballistic. In contrast, the on-set of localization precursor 
effects is pushed towards energies/temperatures lower than the inverse relaxation time 
in the Drude approximation. More specifically, the longer is the range of the impurity 
potential the later localization effects appear. This result may explain why experiments 
have so often failed to detect localization precursor effects in cuprates superconductors. 
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